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Abstract 

In this study, we examine the decoherence of independently coupled qubit systems by using the 
influence functional. We especially concentrated on the effect of qubit flip process. We examine 
the zero-dimensional qubit and one-dimensional qubit systems coupled with the multi-band one- 
dimensional system. The qubit flip process leads to self-excitation and oscillation. We obtained the 
time evolution of von Neumann entropy. We also estimated numerically time evolution of density 
matrix and the entropy. One of the results indicates that the slower random telegraph noise leads 
to the more mixed states. Thus the background charge fluctuation is important for charge qubit. 
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I. INTRODUCTION 



Among various proposals for quantum computations, quantum bits (qubits) in solid state 
materials, such as superconducting Josephson junctions^, and quantum dots^^, have the 
advantage of scalability. Such coherent two level systems constitute qubits and a quantum 
computation can be carried out as a unitary operation applied to many qubit systems. It is 
essential that this quantum coherence is maintained during computation. However, dephas- 
ing is a hard to avoid, due to the interaction between qubit system and the environment. 
The decay of off-diagonal elements of the qubit density matrix shows dephasing. Various 
environments can cause dephasing. The effect of phonons is important in solid systems,-. 
The effect of electromagnetic fluctuation has been extensively studied for Josephson junction 
charge qubits-. The background charge fluctuation is important for charge qubit system J 

The exact expression of Green's functions are obtained for a free fermion and boson system 
and multi-band system.-^ By using these expressions, we consider that qubit systems is cou- 
pled to multi-band one-dimensional bath. Integrating over the multi-band one-dimensional 
system's variables, we obtain the qubit density matrix.— Using this influence functional 
method, we examine the dephasing of the qubit density matrix. In the present study, we 
concentrated on the effect of qubit flip process. The qubit flip process lead to self-excitation 
and oscillation.— We obtained the time evolution of von Neumann entropy. We also present 
the numerical calculation of the density matrix and entropy. One of the results indicates 
that the slower random telegraph noise leads to the more mixed states. Thus the background 
charge fluctuation is important for charge qubit. 

II. HAMILTONIAN 

We examine the Hamiltonian of multi-band fermion or system^ and qubit systems. The 
type of bath and qubit interaction is XY, 
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H mt = ^Tiv^(/+C + /r c j!) (3) 



c ia ,d!£. are bath creation and destruction operator at site i of one-dimensional chain, o is 



mt 

spin index and m is multi-band index. For canonical commutation relation is given by 
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The upper sign is fermion and down sign is boson. = x,y,z) is qubit operator at 

site i,jN is geomagnetic ratio of qubit. Above model, the one-dimensional bath and zero- 
dimensional qubit system interacts with contact interaction. We also examine the one- 
dimensional qubit system where one-dimensional bath interact, the effect of indirect inter- 
action appears. 



III. SINGLE QUBIT SYSTEM 

First, we consider a single qubit and one-dimensional chain system. Thus, we examine the 
effect of direct interaction by influence functional method.— ^^2. The interaction Hamiltonian 
is as below, 

H mt = lN H 2 [A ± (I+c + Iocl)}. (5) 

We integrate about the one dimension bath system. This lead density matrix of qubit 
system, 

rh+(t)=i f z+ j z -(t)=i f z _ 
p(I f z+ ,lL) = I . [dl z+ ][dh-] 



(6) 



I z+ (0)=l[ + ,I z+ (0)=I f z+ 

exp(l(/ g6 [/,_] - I qb [I z+ ]))F[I z+ ,I z _], 



the influence functional is 

F[I Z+ ,I Z _] = [ UScU[dcU[dc +i \[dc. 



(5(c +l (t)- C _ l (t))p( Q (0),c J (0)) 

ex P {^(/[c + ]-/[c_])}, (7) 



where, I g b[I] = Jq^AI z . For following discussion, we choose the action of system, I[c] = 
I [c] + Ii n t[c, S^J, thus the unperturbed is given by, 

/o[c] = (A^r r< p dhdh 

4 Jo Jo 

{cl{ti)A+ 0p - 1 (t 1 ,t,t 2 ,t)c i {t 2 ) 

+ c l (t 1 )A 00 + p - 1 (t 1 ,i,t 2 ,i)c\(t 2 )), (8) 

where A 00p (ti,i,t 2 , j) is the propagator of environmental system which is defined on closed 
time path and has four components 

A 00p (i,ti,j,t 2 ) 

A oo + (^ h) A^(i, ti,j, h) 
Ao^(i,t h j,t 2 ) A^(i,tx,j,t 2 ) r 

(9) 

Introducing the incoming interaction picture for the environment system, we can easily verify 
that Eq. ([7j) turns out, 14 

F[W_] = exp[- / 7jV ^ ±)2 f [ h dhdt 2 
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(I ++ (t 1 )A++(i,t u i,t 2 )I ++ (t 2 ) 
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- I++(t 1 )A+-(i,ti,i,h)I+-fo) 
~ I + 4t 1 )A^ + {i,t 1 ,i,t 2 )I ++ {t 2 )) 

(I„ + (t 1 )A++(i,t 1 ,i,t 2 )I„ + (t 2 ) 
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(10) 

For the convenience we change of the coordinate, 

Vz = {h+ + i*-)/2, £ 2 = (4+ - I z -)/2. (11) 

The 77 and £ are called sojourn and blip. In terms of this variable, density matrix is described 
by as follows 

p(Vz = 1) = I T><T \,p{r)z = -i) = I i><i |, 



p(6 = i) = |TXI,p(6 = -i) = U><TI- 

Other variables are defined by 

r? + = (J ++ + J + _)/2, = (/++ - J+_)/2, (12) 

77_ = (J_ + + J__)/2, = (/-+ - /-)/2. (13) 
In terms of above variables the elements of density matrix are expressed as 
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Therefore, the below equations are hold for these new variables, 

Vz = £+ = 6* = V+ = -V-- (14) 
In this case, the influence functional is expressed as 

FM = eM-i^ 1 f £ dt ldh 
{e(ti)(2G A (t 1 ,z,t 2 ,z))r ? (t 2 ) 
+ r ? (t 1 )(2G R (t 1 ,^,t 2 , ? ))e(t 2 ) 

- 77(ti)(2G fl (ti,i,t 2 ,i))e(* 2 )}], (15) 

where G R (ti, i, t 2 , j), G A (t 1 ,i,t 2 , j) and G K (t\, i, t 2 , j) are retarded Green's function, ad- 
vanced Green's function and Keldysh Green's function. For above integral equation, we 
slice the time and take a difference, we get the differential equation for density matrix is 
given by, 

d(>b(t) ijTT / ,\ i i 



^[H qb ,p h (t)]-^ N HA ± ) 2 <**i 

Pb(ti), (16) 



2G R (i,t,Mi) 
2G A (M,Mi), -2G K {i,t,iM) 



where pb(t) is 



= *w = i, *m = -i 
t(o = i, t(o = -i 



Next we choose the representation of density matrix for the spin diagonal case, 



dp s {t) 
dt 



= ~[H qb , Ps {t)\ - -AiNfrA^) 2 [ dh 



( 



—G K (i,t,i,ti) \* 
G R (i,t,i,tx) + G A (i,t,i,ti) 
i(G A (i, t, i, ti) — G R (i, t, i, ti)) 



\ 
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(18) 



where the density matrix, p s (t) is represented as 



Ps(t) 



/|T(t) ><T(*)I + U(*)X 

\Ut) ><i(t)\, 

11(0 ><T (01, 
VlT(t) ><T (01-11(0 ><H0i/ 



(19) 



The trace of density matrix decreases with time. This event is decoherence. For spin flip 
process, another diagonal element increases with time, this represents self-excitation. One 
of off-diagonal element shows decoherence. Another off-diagonal element shows oscillation 
that modulates the signal. 

IV. ONE-DIMENSIONAL QUBIT SYSTEMS 

Next, we examine one- dimensional qubit systems by using influence functional.— In this 
case, the effect of indirect interaction appears. The interaction Hamiltonian is as below, 



We integrate about spin-less bath system and the qubit systems except the i-th site qubit. 
This lead density matrix of qubit system, 



H mt = lN h 2 J2[A±(^c l + I-cl)}. 



(20) 




exp(-(/g 6 [/ 2i _] - I qb [I zi+ ]))F[I zi+ ,L 



(21) 
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the influence functional is 



F[I zi+ ,I zj _] = J IU[dc^][dcU[dc+i][dc-i] 

5(c +t (t)-c^(t))p(c t (0),c,(0)) 



exp{l(/[ C+ ]-7[c_])}. 



(22) 



where, I q b[I] = J* * hAI iz . For following discussion, we choose the action of system, 7[c] 
7 [c] + I int [I Z i, c], thus the unperturbed is given by, 

(A ±lN h) 2 



4 

(4(t 1 )A+ -- 1 (t 1 ,t,t 2 ,j)c J (t 2 ) 
00 + p - 1 (t 1 ,i,t 2 ,j)c]{ 

(23) 



+ CiitjAak \t 1 ,i,t 2 ,j)c]{t 2 )), 



Introducing the incoming interaction picture for the environment system, the equation turns 
out, 

F[I zi+ ,I Z i-\ = exp[-2 / / dtidt 2 

4 Jo Jo 

(I i++ (t 1 )A+ + (i,t 1 ,j,t 2 )I j++ (t 2 ) 
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- J i+ _(* 1 )AM + (Ml,j,t 2 )Wt2)) 

(7 i _ + (t 1 )A+ + (i,t 1 ,j,t 2 )7 j _ + (t 2 ) 
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- 7 i __(t 1 )A -+(i,t 1 ,i,t 2 )7 j _ + (t 2 ))]. 

(24) 

For the convenience we change of the coordinate, 

Vzi = (/«+ + /«-)/2, e« = " J«-)/2. (25) 
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The i] and £ are called sojourn and blip. In terms of this variable, density matrix is described 
by as follows 

pijlzi = 1) = | U><U \,p(Vzi = -1) = | U><U |, 
P (£ zi = 1) = | U><U \,p{Zzi = -1) = I U><h I- 

Other variables are defined by 

Vi+ = + / i+ _)/2, = (J i++ - J i+ _)/2, (26) 

77,_ = (/*_+ + /i__)/2, 6- = - A— )/2. (27) 
In terms of above variables the elements of density matrix are expressed as 
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Therefore, the below equations are hold for these new variables, 

Vzi = &+ = -&-,€zi = Vi+ = -Vi— ( 28 ) 
In this case, the influence function is expressed as 

F[Vh€i\ = exp[-i / / dtidt 2 
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{U(t 1 )(2G A (t 1 ,t,t 2 ,j)) Vz:j (t 2 ) 
+ Vzi(t 1 )(2G R (t 1 ,t,t 2 ,j))^ J {t 2 ) 
- U(ti)(2G R (t u i,t 2 J))^ zj (t 2 )}], 



(29) 



where G R (ti,i,t 2 , j), G A (ti,i,t 2 ,j) and G K (ti, i, t 2 , j) are retarded Green's function, ad- 
vanced Green's function and Keldysh Green's function. For above integral equation, we 
slice the time and take difference, we get the differential equation for density matrix is given 
by, 



dpb(i,j,t) i i 2 

J t = -p.Hqb,Pb(t)] ~ -{lNhA L ) 



J k 







2G R (i,t,k,h) 



2G A (i,t,k,t 1 ) -2G K (i,t,k,ti] 



Ph(h,k,j). 



Next we choose the representation of density matrix for the spin diagonal case, 

( -G K {i,t,k,t 1 ) \ 



G R (t,t,k,t 1 ) + G A (i,t } k,t 1 ) 
-i(G A (i,t,k,t l ) - G R (i,t,k,t l )) 
G K (i,t,k,t 1 ) 



Ps(k,j,h 



The Green's functions of free bath system are given by^ 

w — \oj(k) — n) + 
1 



G A (k,w 

ik 



w — (uj(k) — fi) — iO' 
G k (k, w) = -2iti{l^2n k )5{w- (w(Jfc) -//)), 
1 



(30) 



(31) 



(32) 



e k B T ±1 

where up sign is fermion and down sign is boson, \i is chemical potential. The trace of density 
matrix decreases with time. For qubit flip process, another diagonal element increases with 
time, this represents self-excitation. One of off-diagonal element shows decoherence. Another 
off-diagonal element shows oscillation that modulates the signal. 



V. MULTI-BAND ONE-DIMENSIONAL SYSTEM 



Next, we examine multi-band qubit systems by using influence functional.— The interac- 
tion Hamiltonian is as below, 



T12 



H int = lN n 2 [M#<%, + J r4 m )]- 



(33) 
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We integrate about the qubit system except the i-th site qubit. m is multi-band index. This 
lead density matrix of qubit system, 



rhi + {t)=i f zi+ ,hi-(t)=i f zi - 

P{lL + ,lL) = / , Udl zi+ ][dl z 

Ji^ ,(o)=iL , .i,i.L(o)=iL 



exp(^(/ gb [J 2i „] - I qb [I zi+ ]))F[I zi+ , I zi _\. (34) 



The influence functional is 

F[I zi+ ,I zj _] = J U^ m [dc^ a ][dc%}[dc^ a ][dc™ ja ] 
^(^, CT (t)-c m li(7 (t))p(c-(0),c-(0)) 
exp{^(/[<J -/[<.])}• 



(35) 



where, I q b[I] = jlhAI iz . For following discussion, we choose the action of system, I[c a ] 
Io[ c a] + Iint[Izi, c a], thus the unperturbed is given by, 

{A LlN hf 



4 

(4 m (*i) A+;-^ m ,(t 1 , i, t 2 ,j)c£(t 2 ) 



'OOp^mcr'm'V" 1 ' "tiJJ^jtr' 

+ C(ii)A - ;;Lw(^i^^2,j>^'(t 2 )), 



(36) 



Introducing the incoming interaction picture for the environment system, the influence func- 
tion is expressed as 

4 Jo 
{77i(*i)G fl (ti, Z, (7, m, t 2 , j, <r, m)^(t 2 ) 

+ Z, (7, TO, t 2 , j, (7, to)to^ 2 ) 

- i] i (t 1 )G K '(*!, i, a, to, t 2 , j, a, to)to(£ 2 )} 
+ {?7i(*i)G fJJ (*i, i, 0", ™, <2, J, o", ^)O(^) 

+ ^(ti)G A (ti, i, (7, TO, f 2 , j, (7, TO)TO(i 2 ) 

- rj i (t 1 )G K : (ti, i, a, to, t 2 ,j, a, m)7fc(* 2 )}] 

\2 /■* /-ti 



expf— - — / / dt±dt2 

* Jo Jo ±! 
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{&(t 1 )(2G A (t 1 , i, a, m, t 2 ,j, a, m))^(t 2 ) 



+ 7] i (t 1 )(2G R (t 1 , i, a, m, t 2 ,j, a, m))^(h 



- Uh) (2G R (t u i, a, m, t 2 , j, a, m))^(t 2 )}), 



(37) 



where G R (ti,i, a, m, t 2 ,j, a, m), G A (t\, i, a, m, t 2 ,j, a, m) and G K (ti, i, a, m, t 2 ,j, a, m) are re- 
tarded Green's function, advanced Green's function and Keldysh Green's function. For above 
integral equation, we slice the time and take difference, we get the differential equation for 
density matrix for the spin diagonal case, 

dPs{iJjt) = ~[H qb , Ps (t)} - {{inHa^y f dtl £ 

k,a,m 



dt 



—G (i, t, a, m, k, ti, a, m) 
G R (i, t, a, m, k, t\, a, m) + G A (i, t, a, m, k, ti, a, m) 
—i(G A (i, t, a, m, k, ti, a, m) — G R (i, t, a, m, k, ti, a, m)) 
V G K (i,t,a,m,k,ti,(T,m) 
Ps(k,j,h). 



\ 



J 



(38) 



Even if multi-band system, we only need diagonal band Green's functions. The Green's 
functions are obtained.- The trace of density matrix decreases with time. For qubit flip 
process, another diagonal element increases with time, this represents self-excitation. One 
of off-diagonal element shows decoherence. Another off-diagonal element shows oscillation 
that modulates the signal. 

The von Neumann entropy is given by following equation. 

d „ 



S = — Trplogp = — lim — — p" 

n-*l on 

The time evolution of the entropy is obtained as below, 



(39) 



dS_ 
~dt 



Tr[p s (j,Mi) 

(l + log(p s (j,Mi))]- 



(40) 
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FIG. 1: Time dependence of trace for {-y N hA ± J c ) 2 =0.01, plt><TI+ll><il(0) = 1 One axis is time 
and another axis is time constant. 

VI. NUMERICAL RESULTS 

Here we present the numerical calculation of the density matrix and entropy. 

When environment is random telegraph noise system, the analytical expression is G K (t) = 
-ill 

iJce T o . The results show with increasing time constant the trace decreases and another 
diagonal element increases. This indicates that the slower random telegraph noise leads 
to the more mixed state. Thus the background charge fluctuation is important for charge 
qubit. 7 

Next, we examine the 0-dimensional free fermion bath. The density matrix under rotating 
wave approximation, four qubit density matrix are given by Fig. 3, 4, 5 and 6. These results 
show analytical described behavior. The von Neumann entropy also calculated as Fig 7. 
The entropy increases with time increases. This quantity is symmetric about origin of time 
and oscillates with changing the chemical potential, because the at t=0 time qubit and 
environment does not entangle. It should be noted that present numerical estimations are 
rotating wave approximation. Thus this numerical study is correct only short time regime. 
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FIG. 2: Time dependence of another diagonal element for (7tv/l4j_ Jc) 2 =0.01, plt><TI llxJ-l(o) = 
1. 

VII. CONCLUSION 

In summary, we had examined decoherence of qubit systems coupled with one-dimensional 
bath. The examined decoherence is the case of single qubit system and one- dimensional qubit 
systems. We obtained the differential-integral equation for general initial condition. The 
trace of density matrix decreases with time. Another diagonal element increases with time, 
this represents self-excitation. One of off-diagonal element shows decoherence. Another off- 
diagonal element shows oscillation that modulates the signal. Decoherence without trace 
conservation occurs. We obtained the time evolution of von Neumann entropy. We also 
present the numerical calculation of density matrix and the entropy. For the random tele- 
graph noise, the results indicate that the slower random telegraph noise leads to the more 
mixed states. Thus the background charge fluctuation is important for charge qubit. 
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FIG. 3: Time dependence of trace for ("y n HA ±_) 2 =0.01, pITxTI+llxil^o) = 1 One axis is time and 
another axis is chemical potential. 

1 Y. Nakamura, Yu. A. Pashkin, and J. S. Tsai : Nature 398 (1999), 786. 

2 T. Hayashi, T. Fujisawa, H-D. Cheong, Y-H. Jeong and Y. Hirayama : Phys. Rev. Lett. 91 
(2003), 226804. 

3 T. Tanamoto : Phys. Rev. A 61 (2000), 22305. 

4 D. Loss and D. P. DiVincenzo : Phys. Rev. A 57 (1998), 120. 

5 T. Fujisawa, T. H. Oosterkamp, W. G. van der Wiel, B. W. Broer, R. Aguado, S. Tarucha and 
L. P. Kouwenhoven : Science 282 (2000), 5390. 

6 A. Shnirman, G. Schon, and Z. Hermon : Phys. Rev. Lett. 79 (1997), 2371. 

7 T. Itakura and Y. Tokura : Phys. Rev. B 67 (2003), 195320. 

8 A. M. Zagoskin : Quantum Theory of Many-Body Systems: Techniques and Application 
( Springer- Verlag New York, 1998) 

9 T. Itakura and N. Kawakami : J. Phys. Soc. Jpn. 64 (1995) 2321. 

10 R. P. Feynan and A. R. Hibbs, Quantum Mechanics and Path Integrals (Mc Graw-Hill, New 
York, 1965) 

11 U. Weiss : Quantum Dissipative Systems (World Scientific, Singapore, 1999) 2nd ed. 



14 



10 



FIG. 4: Time dependence of another diagonal element for (■jNhA±) 2 =0.01, pIT^TI U><il(o) 
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FIG. 5: Time dependence of off-diagonal element for (7at7l4_|_) 2 =0.01, ^ITxll • (q) = 1. 
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FIG. 7: Time dependence of von Neumann entropy for (7tv7l4_l) 2 =0.01, /O s (0) = (1, 1; 1, 1). 
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